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Summary. The Geneva evolutionary code has been modified to study the advanced 
stages (Ne, O, Si burnings) of rotating massive stars. Here we present the results of 
four 20 M0 stars at solar metallicity with initial rotational velocities, Wini, of 0, 100, 
200 and 300km/s in order to show the crucial role of rotation in stellar evolution. 
As aheady known, rotation increases meiss loss and core masses [4]. A fast rotating 
20 M0 star has the same central evolution as a non-rotating 26 Mq . Rotation also 
increases strongly net total metal yields. Furthermore, rotation changes the SN type 
so that more SNIb are predicted (see [5] and [6]). Finally, SN1987A-like supernovae 
progenitor colour can be explained in a single rotating star scenario. 



1 Computer model 

The computer model used is the Geneva evolutionary code (see [5]). Con- 
vective stability is determined by the Schwarzschild criterion. The overshoot- 
ing parameter, aover = dovcr/Hp is equal to 0.1 for Hydrogen- and Helium- 
burning cores and equal to afterwards. Modifications have been made to 
study the advanced stages of the evolution of rotating massive stars. Dynam- 
ical shear has been included using Ric = 1/4 [2]. Note that the computer 
model still includes secular shear and meridional circulation. The structure 
equations have been stabilised using Sugimoto's prescription [7]. Furthermore, 
convection is treated as diffusion from the Oxygen (O) burning stage because 
convection is no longer instantaneous. The algorithm developed for rotational 
mixing is used for this purpose. Finally, the nuclear reaction network has been 
extended and contains all the multiple-a elements up to ^^Ni except ^Be. The 
reaction rates are taken from the NACRE compilation or Hauser-Feschbach 
code calculations (ULB, Belgium). 
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2 Evolution 

The early evolutionary stages are presented in [5]. Here we concentrate on 
solar metallicity 20 stars and study the effect of rotation by examining 
four models with initial rotational velocities of 0, 100, 200 and 300 km/s. 
Calculations have been followed until end of central 0-burning for the Vmi — 
100 and 200 km/s models, end of central Si-burning for the v-mi = 300 km/s 
model and end of first shell Si-burning for the non-rotating model. 

2.1 Hertzsprung Russell (HR) diagram 
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Fig. 1. HR diagram for the solar metallicity 20 M© stars with initial rotational 
velocities of 0, 100, 200 and 300 km/s. 

Figure 1 shows the evolutionary tracks of the four different 20 M© stars 
in the HR diagram. The non-rotating model ends up as a red supergiant 
(RSG) like other group models (see e.g. [1] or [3]). However, the rotating 
models show very interesting features. Although the 100 km/s model remains 
a RSG, the 200 km/s model undergoes a blue loop to finish as a yellow- 
red supergiant whereas the 300 km/s model ends up as a blue supergiant 
(BSC). Thus rotation may strongly affect the shock wave travel time through 
the envelope when the star explodes in a supernova event, since this time is 
proportional to the radius of the star (RSG radii are about hundred times 
BSG ones). Moreover, the behaviour of the models with t^ini between 200 and 
300 km/s is reminiscent of the evolution of the progenitor of 1987A indicating 
that rotation may play a role in similar cases. 
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2.2 Central evolution 

The central evolution is best seen in the central temperature, Tc, versus central 
density, pc, diagram (Fig. 2). We can see that rotation makes the cores slightly 
less degenerate (higher Tc and smaller pc). This is explained by the bigger 
core masses. We also see that the "C-bump" due to the convective central 
C-burning fades away when rotation increases (see also Fig. 3). This is again 
a consequence of more massive cores in rotating models which implies higher 
neutrino loss rates and smaller central carbon abundance at the end of He- 
burning phase. 
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Fig. 2. Tc 


versus pc 


diagram. 



3 "Pre-SN" models 

3.1 Mass loss and core masses 

We can see in Fig. 4 that both mass loss and Helium (He) core masses, M^, 
increase with rotation as already known. There is a saturation effect at high 
rotation when the star is left with hardly any Hydrogen (H) envelope. As can 
be seen in [5] , rotation noticeably increases the number of Wolf-Rayet stars 
(WR). Here we see that there is a smooth transition between SN type from 
HP HL Hb (^ lb). We also note that the Wjni = 300 km/s 20 Mq model 
has a bigger He core than the non-rotating 25 model (it would correspond 
to the core of a non-rotating 26 M© model). The Carbon-Oxygen core mass, 
Mcoi increases with rotation in a similar way as Mq,. The Silicium (Si) core 
mass at the end of central 0-burning only slightly increases with rotation. 



Hirschi et al. 



- 





H 3 A 2 2 » E) A 



2 



Fig. 3. Kippenhahn diagrams; mass (M©) versus logio(~ time left until core col- 
lapse) (yr). Left non-rotating 20 M© model. Right Vmi = 300km/s 20 Mq model. 
Coloured zones show convective zones. Letters indicate burning stages. 
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Fig. 4. Final total mass, Mflnai, and Helium core mass, Ma, as a function of Dini. 
Non-rotating 25 M© model masses are also shown for comparison. 



3.2 Abundances profile and net yields 

In Fig. 5 we notice the smoother profiles due to rotational mixing and also 
the very small quantity of remaining H. The "pre-SN" net yields calculated 
at this stage show that rotation increases the total metal yield and ^^O yield. 
Typically, the total metal yield of the Wini — 300km/s model is twice the one 
of the non-rotating model. On the other hand, rotation decreases H-burning 
products yields (notably ^He) as can be seen in Fig. 6. 
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Fig. 5. Abundances profile for main elements at the end of central Si-burning. Left 
non-rotating 20 M© model. Right Vini ~ 300km/s 20 Mq model. 



M= 20 Mg , Z= & aover= 01 







1 1 1 1 


















^He 


II 




1 , , , '^•i 1 



100 200 300 

"^ini [km/s] 



Fig. 6. Net yields of the sum of all metals, Yz, and individual elements as a function 

of Uini. 
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